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Topological superconductors which support Majorana fermions are thought to be 
realized in one-dimensional semiconducting wires coupled to a superconductor [IHS] ■ 
Such excitations are expected to exhibit non-Abelian statistics and can be used to 
realize quantum gates that are topologically protected from local sources of decoher- 
ence [HIS]. Here we report the observation of the fractional a.c. Josephson effect in 
a hybrid semiconductor /superconductor InSb/Nb nanowire junction, a hallmark of 
topological matter. When the junction is irradiated with a radio-frequency /o in the 
absence of an external magnetic field, quantized voltage steps (Shapiro steps) with 
a height AV = /i/o/2e are observed, as is expected for conventional superconductor 
junctions, where the supercurrent is carried by charge-2e Cooper pairs. At high mag- 
netic fields the height of the first Shapiro step is doubled to hfo/e, suggesting that 
the supercurrent is carried by charge-e quasiparticles. This is a unique signature of 
Majorana fermions, elusive particles predicted ca. 80 years ago [6]. 

In 1928 Dirac reconciled quantum mechanics and special relativity in a set of coupled 
equations, which became the cornerstone of quantum mechanics[7]. Its main prediction that 
every elementary particle has a complex conjugate counterpart - an antiparticle - has been 
confirmed by numerous experiments. A decade later Majorana showed that Dirac's equation 
for spin-1/2 particles can be modified to permit real wavefunctionsP, [8]. The complex 
conjugate of a real number is the number itself, which means that such particles are their 
own antiparticles. While the search for Majorana fermions among elementary particles is 
ongoing [9], excitations with similar properties may emerge in electronic systems[l], and are 
predicted to be present in some unconventional states of matter [TOHTS] . 

Ordinary spin-1/2 particles or excitations carry a charge, and thus cannot be their own 
antiparticles. In a superconductor, however, free charges are screened, and charge-less spin- 
1/2 excitations become possible. The BCS theory allows fermionic excitations which are 
a mixture of electron and hole creation operators, 7^ = c| + Cj. This creation operator is 
invariant with respect to charge conjugation, c| ^ Ci. If the energy of an excitation created 
in this way is zero, the excitation will be a Majorana particle. However, such zero-energy 
modes are not permitted in ordinary s-wave superconductors. 

The current work is inspired by the paper of Sau et al. [15] who predicted that Majorana 
fermions can be formed in a coupled semiconductor /superconductor system. Superconduc- 
tivity can be induced in a semiconductor material by the proximity effect. At zero magnetic 
field electronic states are doubly-degenerate and Majorana modes are not supported. In 
semiconductors with strong spin-orbit (SO) interactions the two spin branches are separated 
in momentum space, but SO interactions do not lift the Kramer's degeneracy. However, in 
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a magnetic field B_LBso there is a range of energies wliere double degeneracy is lifted[T6]. 
see schematic in Fig [lj:;. If the Fermi energy Ep is tuned to be within this single-mode range 
of energies, Ez > -w + Ep, (where A is the proximity gap, Ez = qhbB /2 is the Zeeman 
energy, /is is the Bohr magneton, and g is the Lande ^f-factor), the proximity effect from a 
conventional s-wave superconductor induces p-wave pairing in the semiconductor material 
and drives the system into a topological superconducting state which supports Majorana 
particles. Theoretically, it has been predicted that proper conditions for this to occur can be 
realized in 2D [T^[T7] and, most relevant to the current work, in ID systems [li|2j. Moreover, 
multiband nanowires should also support topological superconductivity P^j - EU] . 

What are the experimental signatures of Majorana particles? Majorana particles come in 
pairs, and zero energy Andreev end-modes localized at the ends of a wire can be probed in 
tunneling experiments [ TTl 1^ ^FT\ . Indeed, there are reports of zero bias anomaly observed in 
topological insulator/ superconductor ^23] and semiconductor/ superconductor j24i [25] struc- 
tures. However, conductivity enhancement near zero bias can also be a signature of diverse 
phenomena in mesoscopic physics, such as the Kondo effect in quantum dots [261 EZ] or the 
"0.7 anomaly" in nanowires [281 129] . Fusion of two Majorana modes produces an ordinary 
fermion and, uniquely to Majorana particles, modifies periodicity of the Josephson relation 
from 27r (Cooper pairs) to An (Majorana particles) [1] B] [301-I32] . In the dc Josephson effect, 
fluctuations between filled and empty Majorana modes will mask the 47r periodicity and, 
indeed, we observe only 27r periodicity in a dc SQUID configuration. In the ac Josephson 
effect, however, the 47r periodicity due to Majorana modes should be fully revealed. 

In our experiments Nb/InSb/Nb Josephson junctions (JJs) are fabricated lithographically 
from a shallow InSb quantum well. Superconductivity in InSb is induced by the proximity 
effect from a Nb film placed on top of the InSb nanowire. The self-aligning fabrication 
process which we use is described in the Supplementary Information. A pattern of multiple 
JJs is defined by e-beam lithography, and a 45 nm layer of Nb is deposited by dc sputtering 
on top of the InSb quantum well. An image of a JJ region is shown in Fig. [1^. A weak 
link is formed between two 120 nm-wide and 0.6 /zm-long Nb wires, with gaps ranging from 
20 to 120 nm in different devices. During deposition of Nb, a thin (2-4 nm) layer of Nb is 
extended 70-80 nm outside the pattern, including the space inside the gaps (a brown halo 
around the Nb wire on the AFM image). This layer is used as an etch mask to define a 
nanowire in the underlying semiconductor self-aligned to the Nb wire. After the etching a 
continuous w ^ 290 nm-wide InSb wire is formed under the Nb, as shown schematically in 
Fig. [1]d. The thin Nb layer is not conducting at low temperatures, so that the supercurrent 
is carried by the proximity-induced superconductivity in InSb. 

The InSb wires have rectangular cross section w ^ d and SO interactions are dominated 
by the Dresselhaus term Hd = '^D{k1){kxax — kyay), where {k"^) = (vr/d)^, d = 20 nm is the 
quantum well thickness, 7^) = 760 eV-A^ for InSb, cTj are Pauli matrices, and x and y are 
the principal crystallographic axes. We estimate Eso ~ 1 meV[33]. The wires are oriented 
along the [110] crystallographic direction, and expected direction of the effective spin-orbit 
magnetic field Bgo is perpendicular to the current, as indicated by the green arrow on the 
AFM image. At high fields {Ez > \J /S^ + Ep) Majorana particles are expected to form 
inside the InSb wire close to the ends of the Nb wires. At these fields the supercurrent is 
dominated by the fusion of two Majorana particles across the gap, which amounts to the 
le charge transfer. From the lithographical dimensions we estimate that only a few (1-3) 
one-dimensional subbands should be populated in InSb nanowires, however, we expect the 
density of states in the nanowires to be modified by the strong coupling to Nb and the actual 



3 




FIG. 1. Devices layout, (a) Optical image of a sample with several devices. A single dc SQUID 
device is outlined with a red oval. On the enlarged AFM image a single Josephson junction is 
shown. The light areas are Nb. A light brown halo around Nb is a thin 2-3 nm Nb layer which 
defines the width of the semiconductor wire after wet etching. The direction of the spin-orbit field 
Bsc, is indicated by the green arrow, (b) A schematic view of the device, the dots mark the expected 
positions of Majorana particles inside the InSb nanowire. (c) Energy dispersion in a material with 
spin-orbit interaction in the presence of magnetic field BJ-Bso- 

number of filled subbands may be larger. 

As devices are cooled down, a series of superconducting transitions — Tcs and are 
observed (Fig. [2]). The first transition, ~ 6.4 K, is for wide areas, Tc2 = 5.8 K and 
= 1.9 K are for the 1 and 0.12 /xm-wide Nb wires, and is for the JJs. From = 1.17 
K for JJS we estimate a proximity gap A ^ 180 /leV and a semiconductor-superconducting 
coupling A ~ 2. 6 A [34j. Lithographically our devices consist of two JJs in parallel, and 
we can measure the ratio of the critical currents in the two arms r = Ici/Ic2 by measuring 
current modulation in a dc SQUID configuration. The ratio r = 7.3 for JJ7 and r > 10 for 
JJ8, indicating that conduction is dominated by a single junction. In the following analysis 
we will treat our devices as containing a single JJ. 

As seen in Fig. |2| V{I) characteristic for JJ8 measured at the base temperature of 20 
mK exhibits a clear supercurrent region {V = 0), with an abrupt appearance of a finite 
voltage. A small hysteresis is observed for the return critical current Ir ~ Ic, characteristic 
of a resistively shunted J J in an intermediate damping regime [33]. Indeed, we measure high 
leakage to the substrate, and estimate the shunting resistance Rg ^ 1 kQ. The measured 
resistance in the normal state is i?^ = 650 Q, and the actual normal resistance of the JJS 
Rn ~ 2 — 6 kQ, consistent with the number of ID channels estimated from the size of the 
wire. The product RnIc ~ 1 mV> 2A indicates that JJ8 is in a clean limit (weak link 
length Leff < ^,1, where ^ is the coherence length and I is the mean free path), a proper 
condition for the formation of Majorana particles. 
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FIG. 2. Characterization of Josephson junctions. Left: V{I) characteristic of a Nb/InSb/Nb 
Josephson junction JJ8 (40 nm gap) measured at 20 mK. In the bottom inset the same characteristic 
is plotted as / vs V, where the black line is an extrapolation from high V. At y < 0.5 mV the 
excess current is clearly seen. The top inset shows the temperature dependence for JJ7 (30 nm 
gap) and JJ8. Right: Normalized differential conductance is plotted as a function of voltage across 
JJ8, with Rj^ = 650 il. The enhancement of the current seen at low V is a signature of Andreev 
reflection. The maximum enhancement is plotted as a function of B in the inset. 
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FIG. 3. ac Josephson effect and Shapiro steps. Left: differential resistance dV/dl (in kil) of 
JJ8 is plotted as a function of the rf amplitude Vrj and dc current / for /o = 3 GHz. The data is 
measured with low frequency (17 Hz) ac excitation lac = 2 nA at T = 20 mK and 5 = 0. |1^(-/^)| 
characteristics at V^j = 0.2, 3, 4, 5 &: 6 mV are shown as white lines; their zero is shifted vertically 
and is marked by dashed white lines. The small vertical green bar indicates the scale of 12 /iV. 
Right: V{I) characteristics of the JJ8 in the presence of B||I measured with Vrf between 3 and 6 
mV in 0.6 mV increments. For i? < 2 T, Shapiro steps with a height = hvrf/2e = 6 jjN are 
clearly observed. For i? > 2 T the plateau at 6 //V disappears, and the first step is observed at 12 
/iV. This doubling of the first Shapiro step is a signature of the ac fractional Josephson effect, and 
is a hallmark of a topological superconductivity. 
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Normalized differential conductance, plotted in the right panel of Fig. |2| shows enhance- 
ment at low voltages. This excess current is a signature of Andreev reflection [36| |37] . Most 
important for our measurements is that the excess current, and thus coherent electron trans- 
port, is observed at high in-plane magnetic fields up to 4 T, as shown in the the inset. 

In the presence of a microwave excitation, phase locking between the rf field and the 
Josephson supercurrent gives rise to constant-voltage Shapiro steps in the V{I) characteris- 
tics ai Vn = nhfo/q, where h is Planck's constant, q is the charge of quasiparticles, /o is the 
microwave frequency, and n = 0,±1,±2... ^38j. Shapiro steps for /o = 3 GHz are shown in 
Fig. [3] At i? = we observe steps with the height AV = 6 fiV, consistent with the Cooper 
pair tunneling (g = 2e). AV scales linearly with /o |33J. The evolution of the steps with 
Vrf is best visualized in the dV/dl plots, where steps with < n < 10 are seen at high 
rf powers. A transition from low to high rf power regime is clearly seen in the dV/dl plot 
near Vrf ^ ^ mV and is marked by a green horizontal line. At high rf powers Vrf > 4 mV 
the evolution of the width of the Shapiro steps A/„ follows a Bessel function pattern as a 
function of power, A/„ = A\Jn{2eVrf /hfo)\, where Vrf is the rf amplitude at the junction. 
We can find the rf power attenuation from the fit to the AJq, Vrf = 5 ■ 10~^Vrf for /o = 3 
GHz, Fig. |4| Here Vrf is the rf amplitude at the top of the fridge. Thus, Vrf = 4 mV 
corresponds to Vrf = 20 /xV ~ Vr/2a/2, where Vr ~ 60 /iV, see Fig. [2] 

For Vrf < 4 mV the junction is in the small microwave signal regime [39]. The linear 
response of a JJ has a singularity at w = ±wj, where flj = 2e/hV is the Josephson frequency, 
and the JJ performs a parametric conversion of the external frequency. The JJ8 is in 
the intermediate damping regime and the V{I) characteristic is expected to become non- 
hysteretic in the vicinity of the first step. Indeed, we observe no hysteresis for Vrf > 1.8 
mV. While nonlinear effects can be present at high I and Vrf, we want to stress that the 
first step at the onset of the normal state is due to phase locking between the external and 
the Josephson frequencies, u = ±uj. 

When in-plane magnetic field B||I is applied, Shapiro steps at = 6, 12 and 18 /iV are 
clearly visible at low fields, i? < 2 T. Steps at 12 and 18 /xV remain visible up to -B ~ 3 T, 
while the step at 6 /iV disappears above B ^ 2 T. The disappearance of all steps above 3 T 
is consistent with suppression of the excess current and Andreev reflection at high fields. 

Quantitative comparison of the width of the Shapiro steps A/„(K/) for different B ex- 
tracted from dV/dl data is plotted in Fig. Q For 5 < 2 T steps at 0, 6, 12, 18 and 24 
liV evolve according to Bessel functions with the amplitude A f» 150 nA. For B > 2 T the 
step at 6 /iV vanishes at low Vrf, and re-appears at high Vrf > 10 mV for B = 2.2 T. We 
also observe that the low field rf attenuation does not fit the evolution of the A/q plateau. 
Moreover, a minimum of the 12 fiV plateau at Vrf ~ 7 meV coincides with the minimum of 
the I Ji| Bessel function, suggesting that indeed the 12 /xV plateau became the n = 1 Shapiro 
step. Evolution of plateaus at -B > 2 T is poorly described by Bessel functions, suggesting 
that oscillations with different frequencies may contribute to the width if the plateaus. We 
emphasize that the doubling of the Shapiro step height is a unique signature of a topological 
quantum phase transition. 

Theoretically it has been argued that Josephson current with both 2tt (/csin(0)) and 
47r (Jm sin(0/2)) periodicity should be present in the topological state [3l II0I - H2] . especially 
in a multichannel wires. However, in current-biased junctions odd steps are expected to 
vanish even in the presence of large supercurrents carried by the charge-2e quasiparticles 
Ic ^ Im [13]. In this case the width of even steps is expected to be defined primarily 
by Jc and the coefficient A remains almost unchanged across the transition, as is observed 
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FIG. 4. Evolution of Shapiro steps with rf power. The width of the first five Shapiro steps 
are plotted as a function of the rf field amplitude Vj-j. Lines are Bessel functions A\Jn{f3Vrf)\, with 
/? = 0.84 mV~i (solid) and /3 = 1.04 mV-^ (dashed). 



experimentally, especially if a large number of subbands is occupied in the nanowire due to 
the strong coupling to the superconductor. At high voltages across the junction we expect 
enhanced mixing between gapless and gapped modes, which may explain the prominence of 
the 18 /iV and higher plateaus at -B > 2 T. 
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I. ANALYSIS OF THE PARAMETER SPACE TO OBSERVE FRACTIONAL 

JOSEPHSON EFFECT 

In order to form Majorana fermions in a nanowire, several conditions have to be satisfied. 



The most stringent is lifting of the Kramers degeneracy, Ez > yA^ + Ep, where Ez is 
Zeeman energy Ez = g* hbB, g* = 50 for InSb, fiB is Bohr magnetron, and B is an external 
magnetic field. At the same time we need the proximity gap to be larger than the Josephson 
frequency, A > huj = 6 /xeV for 3 GHz, or A > 2hijjj in the topological phase. The 
proximity gap A depends on the semiconductor-superconductor couphng A, Zeeman energy 
and spin-orbit coupling [34J: 

A = A. ^ \ ^ ^gg^ . (SI) 



The maximum A cannot exceed the superconducting gap in narrow Nb wires, A^ = 290 
fieV for = 2 K, and is maximized for large coupling A ^ A^. Large coupling, though, 
increases electron density in the semiconductor which, in turn, increases Ep, thus requiring 
large fields to lift the Kramers degeneracy. 
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FIG. S5. Left: Fermi energy Ep and energy of spin-orbital coupling Eso (for 20 nm InSb quantum 
well) are plotted as a function of carrier momentum k. The right axis shows magnetic fields which 
correspond to the Zeeman energy of the left axis. Right: expected field dependence of the proximity 
gap (Ed. [Sl| for Eso = 1 meV and A(0) = 180 ^ueV. 



Let us analyze the parameters of the JJ8 device. We measured = 1.17 K, thus 
A = l.TGksT/e = 180 fieV and A = 470 /leV, or A ~ 2.6A. In our sample geometry we 
cannot measure electron density independently, but we can assume that aX B = 2 Tesla 
(observed topological phase transition) Ez = Ep. We plot Ep = {hkY /2m* , E^ = qhbB /2 



and Eso = V^loik^k as a function of k in Fig. S5, Here we use 'Jd = 760 eV-A^, 
(^z) — {'^/dYi effective mass m* = O.OlSmo, effective (^-factor g = 50, and the quantum well 
thickness d = 20 nm. At 5 = 2 T {Ez = 2.8 meV) the condition Ep = Ez translates into 
k ~ 3.5 • 10^ cm~^. Spin-orbit coupling for this k is estimated Eso ~ 1 meV. The expected 
field dependence of the proximity gap A(i?) is plotted on the right plot for Eso = 1 meV. 
We see that at 5 = 2 T A = 30 fieV and the condition A > 2huj is satisfied. 

We observe Shapiro steps up to ~ 3 T, which is consistent with the theoretically expected 
gap to be only 20 fieV at that field. Experimentally we measure a weaker 1^ vs B dependence 



than the one predicted by Eq. (SI), see Section IV 



The maximum width of the InSb wire is 290 nm, and we expect the actual InSb wire 
width to be reduced due to side etching and surface depletion. For an InSb wire with width 
w = 100 — 250 nm the energy separation between the first 2 energy levels Sh"^ /8m*w'^ is 
1.2-8 meV and only < 3 subbands are expected to be occupied for Ep < 3 meV. We note, 
though, that strong coupling to Nb modifies the density of states in the InSb wire and the 
actual number of filled subbands may be larger. Thus, we conclude that experimental 
parameters for the JJ8 device satisfy the requirements for the observation of 
Majorana fermions. 



II. FABRICATION OF JOSEPHSON JUNCTIONS 



The starting material is undoped Ino.6Gao.4Sb/InSb/Ino.6Gao.4Sb (3nm/20nm/3nm) 
quantum well grown by molecular beam epitaxy on a Te-doped (001) GaSb substrate. 
A thick graded In^-Gai-ajSh {x = .17 — .6) buffer and a 120nm Ino.77Alo.33Sb barrier were 




FIG. S6. Optical image of a device and AFM micrographs of a Josephson junction before and after 
wet etching. Yellow dashed line outlines the extent of the thin Nb layer. 



grown between the wafer and the quantum well for strain relaxation and electron confine- 
ment. A pattern of multiple JJs is defined by e-beam lithography, and a 45 nm thick layer of 
Nb is deposited by dc sputtering. Surface oxide is removed by dipping the sample in HF:DI 
(1:10) for 20 seconds prior to Nb deposition. An optical image of a multi-device sample is 
shown in Fig. S6 Weak links are formed between two 120 nm-wide and 0.6 /im-long Nb 



wires, with gaps in 20-120 nm range. 

The key processing step is self alignment of the Nb and InSb wires. We use double layer 
MMA/PMMA photoresist which creates an undercut after e-beam exposure and develop- 
ment. During Nb sputtering at 0.1 mTorr of Argon some Nb is scattered into this undercut 
area and a thin (2-4 nm) layer of Nb is formed, extending 70-80 nm outside the pattern, 
including the space inside the gaps (a brown halo around the wire in the AFM image; see 
also Fig. S7). This layer is used as an etch mask to define a nanowire in the underlying 
semiconductor, which becomes self-aligned to the Nb wire. Etching in H2S04:H202:H20 
(1:8:1000) for 30 sec removes 60 nm of semiconductor and a continuous ~ 170 nm-wide InSb 
wire is formed under the Nb (we expect that the width of the wire is reduced during the 
wet etching step by ~ 60 nm from each side). 



III. PROXIMITY EFFECT 

In order to verify that we observe proximity-induced superconductivity, we fabricated 
several (> 10) test devices on a semi-insulating GaAs substrate with the Nb pattern identical 



to the J J devices, see inset in Fig. [ST) Continuous wires show the expected superconducting 
phase transitions at Td = 7.5 K for wide regions (> 6/im), Tc2 = 7.1 K for 1 /im wide 
connectors and Tcs = 1.9 K for the 80 nm-wide wire. There is a 4 nm thick wetting layer 
around Nb, which can be seen as a halo around wires on the micrographs. For wires with 
small gaps < 100 nm the wetting layer fills the gap. The wetting layer is not attacked by 
the etching solution and serves as an etch mask in the semiconductor wire definition. In 



the device shown in Fig. S7 the gap is ~ 40nm. Yet these devices become insulating when 
cooled to low temperatures. This test experiment allows us to establish that (i) the thin 
Nb layer is not conducting and plays no role in the electrical transport, (ii) the tunneling 
current for gaps > 20 nm is negligible, and (iii) in InSb JJs the current has to flow through 
the InSb layer. Thus the observed superconductivity in InSb JJs is due to the 
proximity effect. 
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FIG. S7. Temperature dependence of a Nb junction and a continuous wire fabricated on a semi- 
insulating GaAs substrate. The wires are fabricated identically to the JJ on InSb substrates. In 
the bottom inset AFM micrographs of the junction and the wire devices are shown (2 /im x 2 ^m 
scan size) . On the top image a zoom of the gap region is shown. A 4 nm-high wetting layer around 
the wire and within the gap is clearly seen and is outlined with a dashed yellow line. The wire 
profile along the white line is shown on the right, and the thin Nb layer is marked as WL. 

IV. TEMPERATURE DEPENDENCE OF JJ 



Proximity-induced superconductivity has a reduced gap A compared to the gap of the 
host superconductor, Eq. (SI). The lowest corresponds to this reduced proximity gap. 

Temperature dependence of samples resistance is shown in Fig. S8 For a continuous 
line (L) we can identify three transition temperatures: Tqi = 6.3 K is for wide regions, 
Tc2 = 5.7 K is for 1 /im-wide wires and Tcs = 1.9 K is for the 150 nm-wide wire. Similar 



results are obtained for Nb on GaAs, Section III Devices with a Josephson junctions ( JJ6-8) 
have the actual superconducting transition < Tc < Tc3 for various devices. Note that in 
a '^He refrigerator, where electrical noise is higher, we do not observe the superconducting 
transition for the JJ8 device down to 250 mK, the actual Tc for this device is 1.17 K, see 
inset. 



V. MAGNETIC FIELD DEPENDENCE OF CRITICAL CURRENT 



The requirement is that the proximity gap A > ksT, hujj in the semiconductor material at 
the magnetic fields needed to satisfy the condition Ez > Ep- In our Nb wires superconduc- 

In the low field region Ic{B) 



tivity survives up to 5 Tesla, see Fig. 



S9 



/,(0)/^l - iB/B,y 

with Bi ^ 2.5 T. Resistivity of crystalline Nb is 152 nfi-m and a perfect L = 1.2 fim, w = 100 
nm and t = AO nm wire should have resistance of ~ 1 fi. Our wires have resistance of 2 — 3 
kfl, which indicates a substantial degree of disorder. Neither A nor are expected to be 
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temperature (K) 

FIG. S8. Temperature dependence of resistance for the devices JJ6 (20 nm gap), JJ7 (30 nm gap), 
JJ8 (40 nm gap) and a continuous 150 nm-wide line (L) is measured in '^He system (main plot) 
and in a dilution refrigerator (inset). The dilution refrigerator is properly shielded. 



L (line) 
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FIG. S9. Magnetic field dependence of differential resistance for device JJ8 and a 150 nm wide 
wire fabricated next to it (L). Black regions are the superconducting states. Yellow line is 
Ic/y/l - {B/BiY, where Bi = 2.5 T is the first critical field. 



substantially affected by the disorder [35] , consistent with the measured Tc ~ 7 K compared 
to the Tc = 9.2 K in pure crystaUine Nb. Disordered Nb is a type-II superconductor, and 
superconducting gap survives to much higher field i?2 ~ 5 T. The thickness of the film t = 40 
nm is of the order of the coherence length in Nb (39 nm), thus flux capturing is possible 
even for the in-plane field. 
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VI. FREQUENCY DEPENDENCE OF SHAPIRO STEPS 



In Fig. SIO we show V{I) traces measured in the presence of rf field with frequencies 

hfo/2e = 4, 



/o = 2, 3 and 4 GHz. The corresponding step heights are, respectively, AV 
6 and 8 /iV. 



f = 4 GHz 
AV= 8 



f = 3 GHz 
AV= 6 |iV 



f=2GHz 
AV= 4 jiV 
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I(nA) 



I(nA) 



FIG. SIO. Shapiro steps for rf frequencies /o = 2, 3 and 4 GHz. 



VII. ANALYSIS OF THE JOSEPHSON JUNCTION 

The wafers used in these experiments have substantial leakage between the Nb film and 
a heavily doped substrate, and the substrate forms a shunting resistance to every device. 
We can estimate the shunting resistance Rg and the normal resistance of the JJ8 as 
follows. Next to the J J we fabricate a 0.12 /xm-wide wire with a measured resistance in 
the normal state of {R~^ + -R^J)"^ = 0.6 kQ. An identical wire on an insulating GaAs 
substrate has a resistance of R^ = 2.5 kQ, thus the shunting resistance for the wire device 
Rg^ ^ 0.8 kQ. The wire device and the JJ devices have similar layouts and we assume that 
the shunting resistances are similar, Rs ~ Rgw (the leakage to the substrate is primarily 
through contact pads of the same size). The measured resistance of the JJ8 in the normal 
state is {R]^^ + Rj^)~^ = 0.65 kf2, thus we estimate Rn ~ 3 kil. The JJ8 device consists of 
two nominally identical JJs in parallel, and the actual Rn for the dominant junction can be 
as high as ~ 6 kfl. Thus, JJ8 has a few (< 5) subbands, which is also consistent with the 
number of subbands estimated from the size quantization, see Section |T} 

For short weak links in the dirty regime {Lgff <C ^ and / <^ -^e//, where / is the mean free 
path, ^ is the coherence length, and Leff is the effective length of the weak link), the product 
/ci?Af[yuV] ~ 320Tc[K] [H]. For the clean regime {Leff ^ the value should be almost 
twice as large, /c-RAr[/iV] ~ 480Tc[K]. From the above estimates the product IcRn ~ 1 mV, 
while Tc = 1.17 K. Thus, JJ8 is in the clean regime, which is favorable for the observation 
of Majorana fermions. 
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FIG. Sll. Left: differential resistance dV/dl of the JJ8 is plotted as a function of the rf amplitude 
Vrf and dc current /. The rf frequency is /o = 3 GHz. The data is measured with low frequency 
(17 Hz) ac excitation lac = 2 nA at T = 20 mK and B = 0. 1^(1)1 characteristics at Vrf = 0.2, 
3, 4, 5 & 6 mV are shown as white lines, their zero is shifted vertically and is marked by dashed 
white lines. The short vertical bar is a 12 //V scale. Right: n-th step width AIn is extracted from 
the dV/dl data and plotted as a function of Vrf for n = — 4. Solid lines are Bessel functions 
A\Jn{pVrf)\ with A = 150 nA and /3 = 0.89 mV"^ A red dashed line forn = 1 is oc Vrf /I. 



VIII. ANALYSIS OF SHAPIRO STEPS AS A FUNCTION OF RF POWER 



In Fig. Sll we plot differential resistance dV/dI{I,Vrf) ai B = for the JJ8 device. 
Shapiro steps with AV^ = 6 /iV are clearly seen at K/ > 2 mV and correspond to the black 
regions on the dV/ dl contour plot. At high Vrf > 4 mV evolution of the width of the Shapiro 
steps AIn follows the Bessel function pattern as a function of power, A/„ = A\Jn{2eVrf / hfo) \ , 
where Vrf is the rf ampUtude at the junction and /o is the rf frequency. We can estimate the 
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rf power attenuation from the fitting of AIq with Jq above Vrf > 4: mV, Vrf = 5.4 ■ 10~^K/ 
at /o = 3 GHz and Vrf = 1.7 ■ 10"^^./ at /o = 4 GHz. Here Vrj is the rf amphtude at the 
top of the fridge. 

In the dV/ dl{l, Vj-j) plot we can identify two regions with different A/o(K-/) dependencies, 
separated schematicaUy by a green hne at V^p^ ~ 3.7 mV. For the JJ8 device Vc = IrRn ~ 
65 fieV ^ hfo/2e = 6/ieV, see Fig. 2 of the manuscript, and for Vrj < V^p^ the junction is 
in the small microwave signal regime [39]. The linear response of a J J has a singularity at 
Lo = ±ujj, and the first Shapiro step appears due to phase locking of the external frequency 
and the Josephson oscillations. The JJ8 is in the intermediate damping (resistively shunted) 
regime, and the V{I) characteristic is expected to become non-hysteretic in the vicinity of 
the first step. Indeed, we observe no hysteresis for Vrf > 1.8 mV. While nonlinear effects 
can be present at high / and Vrf we want to stress that the first step at the 
onset of the normal state is due to phase locking between the external and the 
Josephson frequencies, u = 

The width of the first step in the low rf power regime is expected to be A/i = 2/c/^/7; 
see dashed red line on A/i(V^j) plot. At high power the width of the n-th step is expected 
to follow the n-th Bessel function. Indeed the step widths follow the expected pattern, but 
some deviations are expected due the sample being in a crossover regime between the low 
and the high rf power regimes. 



IX. MAGNETIC FIELD DEPENDENCE OF THE SHAPIRO STEPS 



A detailed field dependence of step evolution is given in Fig. S12 and the analysis of step 



width as a function of rf power is discussed in Section VIII 
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FIG. S12. Evolution of Shapiro steps in JJ8 with magnetic field. The data is measured with lac = 2 
nA at T = 20 mK. The steps at 6 //V and 12 /xV are outlined with red and green, respectively. We 
use V[I) data to identify the step height. Note that the 6 /xV step disappears for i? > 2 Tesla. 



